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Abstract 



This paper continues a study on Choptuik scaling in gravitational collapse 
of a massless complex scalar field at the threshold for black hole formation. 
We perform a linear perturbation analysis of the previously derived complex 
critical solution, and calculate the critical exponent for black hole mass, 7 ~ 
0.387111 ± 0.000003. We also show that this critical solution is unstable via 
a growing oscillatory mode. 
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I. INTRODUCTION 



The recent discovery by Choptuik and others of new and unexpected critical behav- 
ior in gravitational collapse has sparked efforts, both numerical and analytic, to understand 
these phenomena. One would especially like to understand the nature of the black hole 
scaling relation, and how the observed echoing arises out of the Einstein equations. That 
these phenomena appear for a variety of different matter fields suggests that they may re- 
flect universal properties of the Einstein equations. Choptuik originally conjectured that 
the exponent, 7, in the threshold scaling relation for black holes mass, 

1 \ f 0, p < p* 

M BH(p)oc| (p _ pT (7Rj037); p>p * (1) 

was universal for collapse of a scalar field, and presented considerable evidence. Here p 
is some parameter which characterizes the strength of the initial conditions, and p* is the 
threshold value, i.e., the value for the critical solution. Work by Abrahams and Evans on 
vacuum collapse in axial symmetry, and Evans and Coleman on collapse of a radiation 
fluid, suggested a broader universality or near-universality of 7 (see also ||) among different 
matter models. All these results were obtained by numerical relativity, so that accuracy of 
7 was at best 10~ 3 . However, Maison || showed that fluid collapse models with an equation 
of state given by p = kp have a critical exponent 7 that depends strongly on the parameter 
k in the range < k < 0.88. (Evans and Coleman had examined only the case for radiation 
fluid, k = 1/3.) Maison used an accurate perturbation theory method, as suggested by 
Evans and Coleman, and carried out by Koike, Hara & Adachi Perturbation theory 
affords both a calculation of 7, and a test of stability. 

Just at the critical value p = p* — at the "black hole threshold' — a critical solution 
of the field equation appears, a "choptuon". A choptuon is a delicately poised dynamical 
state of the fields, a collapsing, shrinking, radiating ball of field energy. Its self-gravity 
gives rise to spacetime curvature that partially traps field energy, but only partially: Energy 
continually leaks outward and becomes outgoing radiation. By careful tuning of p, the rate 
of collapse exactly balances the radiation rate, so that the collapse becomes self-similar, 
and the choptuon shrinks all the way down to the Planck scale. This situation is of course 
intrinsically unstable: If self-gravity is a little too weak, the choptuon will soon dissipate 
completely; if a little too strong, it will soon form an event horizon and undergo terminal 
gravitational collapse into a black hole much larger than the Planck mass. Construction of 
the choptuon implies stabilization, by careful tuning of a parameter p, against this threshold 
instability toward either subcritical dissipation or supercritical black hole formation. We will 
call this instability the "black hole instability" for short, bearing in mind that its other side 
is dissipation. The critical exponent 7 measures the strength of the black hole instability. 

Then the question of stability of the choptuon is this: Does the choptuon have any other 
instabilities, aside from the inevitable black hole instability? The previous evidence is that 
choptuons are in fact stable in this sense 

In this paper we return to the gravitational collapse of a massless complex scalar field in 
spherical symmetry, as a good theoretical laboratory in which to address these issues. In || 
- hereinafter called paper I — we derived a new continuously self-similar choptuon at the 
threshold of black hole formation, and constructed the complete spacetime. The discrete 
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echoing seen in the original choptuon |[] reappeared as continuous phase oscillations of the 
complex scalar field. However, the choptuon by itself does not yield the black hole scaling 
relation; one needs to perturb the choptuon. In addition, perturbation theory can probe 
important questions about the stability of choptuons, and the nature of the relative attrac- 
tor in the space of all initial conditions. This paper considers only spherically symmetric 
perturbations, the relevant ones for Eq. ([!]), and the most likely ones to be unstable. It 
would also be interesting to carry out a nonspherical perturbation analysis of any of the 
known choptuons, to check Choptuik's "no hair" conjecture 

The properties of the critical solution are described in paper I, and will not be repeated 
here. Since there is confusion in the literature on this point, we emphasize again that the 
self-similar hypothesis — whether discrete [p]] or continuous fHH^ - is not unphysical for 
gravitational collapse in asymptotically flat spacetime. No self-similar spacetime is asymp- 
totically flat when extended to infinity (except flat spacetime); however this is no problem, 
because the self-similar spacetime can be matched onto an asymptotically flat spacetime at 
some large radius, and within this radius the self-similar approximation is a good one; in 
fact, it becomes increasingly good and natural as the collapse proceeds to small scales at 
threshold 

Section II reviews the calculation of critical exponents, while Sects. Ill and IV give the 
perturbation equations and their boundary conditions. Section V deals with gauge modes, 
while the numerical methods are outlined in Sect. VI. 

Our main results are given in Sect. VII. Upon perturbing the critical solution for the 
complex scalar field, we find two growing modes. The first mode is real; as expected it 
corresponds to formation of black holes off threshold, and it yields a value 7 ~ 0.387111 ± 
0.000003. The second growing mode — actually a conjugate pair of modes — is oscillatory 
and must correspond to an instability of the critical solution on threshold. 

Finally Sect. VIII discusses significance of the results. We do not attempt to follow 
the growth of the second growing mode outside the range of validity of linear perturbation 
theory, but plausibly it would develop into the Choptuik discretely self-similar solution, the 
original choptuon 

II. CRITICAL EXPONENTS AND PERTURBATION THEORY 

We briefly review the perturbation theory method |3|-|5| and then apply it to our model. 
Begin with a continuously self-similar choptuon. Perturb it. Quite generally, perturbed 
self-similar solutions to partial differential equations show power-law behavior of modesQ in 
time t, tantamount to exponential behavior in logarithmic time r = log(— t). Consider a 
perturbation mode with power-law behavior e(— t)~ K , where e (oc p — p*) is a small constant, 
and t is a suitable time coordinate, in this case proper time along the time axis r = before 
collapse at t = 0. Assume it develops into a small black hole. The actual formation of 
the black hole is outside the scope of perturbation theory, but we can still determine the 



In degenerate cases there may be logarithmic corrections to power-law behavior, of the form 
e(— t)~ K log n (— t), but these turn out to be absent in the problem at hand. 
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black hole mass by perturbation theory, as follows. By scale invariance of the choptuon, the 
formation of black holes belonging to different values of e will be strictly homothetic to each 
other: i.e., the full spacetimes will be related to each other by exact scale transformations. 
Determine the mass M of the black hole as follows: Pick a fiducial amplitude of the pertur- 
bation mode, say small 5, which is still within the range of validity of perturbation theory. 
Then imagine evolving further by an exact calculation — but the exact calculation need not 
be done, because it will just scale. The black hole mass Mbh will be proportional to the 
value t\ at which 

<~ti)- K = S (2) 

or 

M BH = const(-ti) = const'e 1/K . (3) 

The critical exponent, Eq. (p]), can then be read off as 

7 = V«. (4) 

To obtain the critical exponent 7, it thus suffices to find the appropriate unstable mode in 
perturbation theory, and its characteristic growth rate k. 



III. EQUATIONS OF MOTION 

To perturb about a self-similar critical solution (choptuon), it is convenient to adopt 
coordinates r = log(— t), z = —r/t where r is the usual areal coordinate that measures the 
area of 2-spheres in spherical symmetry, and t is an orthogonal time coordinate, chosen to 
agree with proper time along the time axis (t < 0,r = 0). With these coordinates, the 
time-dependent spherically symmetric metric can be written, 

ds 2 = e 2r ((1 + u) [-{b 2 - z 2 )dr 2 + 2zdrdz + dz 2 ] + z 2 dQ 2 ) . (5) 

To enforce regularity along the time axis z = at the center of spherical symmetry, the 
metric functions w(r, z) and b(r, z) must obey boundary conditions 

b(r, 0) = 1, w(r, 0) = (6) 

The perturbed fields are then defined as 

b(r,z) ^b (z)+eb 1 (r,z) (7a) 
u(t, z) « Uq{z) + eui(r, z) (7b) 
0(r, z) w e^faiz) + ee^<Mr, z) (7c) 

where denotes the 0th order critical solution, 1 denotes the 1st order perturbation, u = 
1.9154446 is the (unique) eigenvalue of the unperturbed equations (paper I); and where 
e > is an infinitesimal constant, a measure of how far away the solution is from the critical 
solution in the space of initial conditions. Using Choptuik's terminology, we consider the 
supercritical regime for infinitesimal 
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e oc p — p*. 



(8) 



Matter consists of a free, massless, 2 complex scalar field <fi that obeys the scalar wave 
equation 



□0 = (9) 



while the Einstein equations are 



Riw ~ ^9nuR = 8vrT^ (10a) 



SttRo ( V M 0V,0 - i^V0V p ) : (101)) 



paper I gives these equations in detail. 

We now perturb these Einstein-scalar equations through 1st order in e, to obtain a set 
of linear partial differential equations for the perturbed fields b\, u\, (pi, in the independent 
variables r, z. Following the standard approach, we Fourier transform the 1st order fields 
with respect to the ignorable coordinate r = log(— t), 



Ui(a,z) = j e lul ~U\{t ; z)dr ; (11a) 
h((T,z) = J e i<TT b x {r,z)dr, (lib) 
^(a,z) = [ e^Mr,z)dr; (11c) 



throughout, "will denote such a Fourier transform. The transform coordinate a is in general 
complex. The 1st order field equations now become ordinary differential equations (ODE's) 
in z, and under appropriate boundary conditions, become an eigenvalue problem for a. So- 
lutions of the eigenvalue problem are then normal modes of the critical solution. Eigenvalues 
in the lower half plane Imcr < belong to unstable (growing) normal modes. Eigenvalues 
in the upper half a plane — which we will not explore in this paper — would correspond 
to quasi-normal (dying) modes of the critical solution. The eigenvalue a is related through 
Eq. (fD to the critical exponent by 7 = — 1/Imcr. 

Thereby, the Einstein-scalar equations reduce to the following set of ODE's in z. Define 
some auxiliary fields in 0th order 

q{z) = M*)' (12) 

p( z ) = t-(«*>0o(*) - zM z Y) ( 13 ) 

and in 1st order 



2 In view of previous results Q it is unlikely that a nonvanishing rest mass, or more generally a 
nonlinear self-coupling of the form □(/> = V'{(j)) would make any difference below. However, other 
self-couplings such as in the axion-dilaton system |7| may well make a difference. 
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r(a,z)=fa(<r,z) (14) 

s((T, z ) = -f ~ <7 )<?i ~ #'i)) ( 15 ) 

i2(<r,z) = r(-<7*,z)* (16) 

S(<7,*)=s(-<7*,z)*. (17) 

where * denotes complex conjugate. The field equations become (in notation which closely 
follows that of paper I), in Oth order, 

b = , (18a) 

z 

u = Attz 2 (\q\ 2 + \p\ 2 + 2jRe(qp*)^ , (18b) 

u' = -Anb (l + u )(qp* +pq*), (18c) 

, (uq + 2 z(3+\ bo/3- , 

, b p+ fu z0-\ 

p=-£r«-[T + -£r) p > (18e) 

and in 1st order 

, (2zb b 1 uJl} ( 2b x b\ iab x + u 1 b \ 
r = {-^ (3+ -^A) q+ {-^- + A ) p - 

(^ + Z -^ K ^)r + b fW--ia)s, (19a) 

, / 26i^ 2 u 1 b \ [2b 1 b z zu x . zb x \ 

— {/3 + -ia)r-\— + Is, (19b) 



J%±2 +£ (^)j B+ | (/3 ,_^ (19c) 
, / 26iz 2 (2b x b Q z zu x . zb x \ 

£ w _ Iff)fl _^ + i(£z«i) s , (1M) 

&i = -(ui6o + «o&i), (19e) 

&i = ~ Ml(1 + M0) + ^ + 4tt(1 + Mo )^ f-^b| 2 + P S + P *s + q*r + qr) . (19f) 



Here 
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, d . 
= s , (20a) 

p± = itu + u ±l, (20b) 
A = 6g - (20c) 

The equations for R and S are identical to those for r and s under the conjugacy 
(q,p,r, s, f3±) — ► (q*,p*, R, S, @±). In the 0th order problem, the Bianchi identities allow 
us to determine % as an algebraic expression, Eq. (|T8|b). Similarly U\(z) can be written as 
a 1st order algebraic expression involving other 1st order variables: 

ui = 47TZ ^° + X) f-^bl 2 + z(qR + Q*r +PS + p*s) + b (pR + p*r + qS + q*s)) . (21) 



IV. BOUNDARY CONDITIONS 

Now we specify boundary conditions for these equations. On the axis of spherical sym- 
metry, z — 0, solutions must be regular. This demand leads to boundary conditions in 0th 
order 

6o(0) = 1, (22a) 

«o(0) = 0, (22b) 

g(0) = 0, (22c) 

p(0) = free real const; (22d) 

(22e) 

and in 1st order 

Si(0) = 0, (23a) 

ui(0) = 0, (23b) 

r(0) = 0, (23c) 

R(Q) = 0, (23d) 

s(0) = free complex const, (23e) 

5(0) = free complex const. (23f) 

The global U(l) phase symmetry present in the unperturbed equations has allowed us to 
set p(0) to be real. No such rotation can be performed for the time-dependent quantity 0i, 
so s(0) and 5(0) are in general complex. 

The field equations have a regular singular point at the point z = z 2 where A(= 6q — z 2 ) 
vanishes. As discussed in paper I, this locus corresponds to a similarity horizon in spacetime, 
the past light cone of the spacetime singularity at the spacetime origin (t,r) = (0,0). All 
fields must be regular on this horizon, because it lies within the domain of dependence of 
the initial hypersurface. This means in particular that the terms proportional to 1/A in 
Eqs. (|19|ab) must cancel at z 2? giving a linear relation between r(z 2 ), s^), ^1(^2) and ^1(^2); 
and similarly for R{zq) and Sfa) from Eqs. (|l9|cd) : 
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(uo - 1) ((J3+ - ia)r - - ia)s) = q + {iuj + 2)u x + b x u' Q 

( vJh(JL - ia) iah ; , \ ; , , 

p I 1 — - — V hu + lUJUx I + »i(p + g - p^p ) 

(uo - 1) ((/?; - ia)R - (/T - ia)S) = q* + ( 

V 1 1- hu - lUJUx I + bi(p+q ~ P-P )■ 

These equations may conveniently be used to solve for s(^) and Siz^) in terms of other 
boundary conditions as long as a lies in the lower half plane; if one wished to explore the 
upper half plane this procedure would fail somewhere due to the vanishing of the coefficients 
of s and S dX z 2 . The resulting boundary conditions at z = z 2 are as follows, in Oth order: 

bo(z 2 ) = z 2 = free real const, (25a) 

p{z2) — free complex const; (25b) 

(25c) 

and in 1st order: 

i>i(z 2 ) = free complex const, (26a) 

U\(z 2 ) = free complex const, (26b) 

r(z 2 ) = free complex const, (26c) 

R(z 2 ) = free complex const. (26d) 

The Oth order equations, Eqs. (|T8|), amount to five real ODE's in five real fields (b,p, q), 
where the complex fields p and q each count as 2 real fields. Correspondingly there are 5 free 
real constants: (u,p(0), z 2 ,p(z 2 )), where the complex p(z 2 ) counts twice.Q Thus we have a 
well posed eigenvalue problem on a finite domain < z < z 2 , and we may expect a discrete 
spectrum. 

The 1st order equations, Eqs. ([19]) amount to 6 complex linear ODE'sQ in the 6 complex 
fields (bi,ui,r,s,R,S). At our disposal are (a, s(0), S(0),bi(z 2 ),ui(z 2 ),r(z 2 ), R(z 2 )), 7 free 
complex constants. Do we have one too many free complex constants? No, because solutions 
to the equations must come in 1-complex-dimensional linear families, due to linearity; that 
is, each solution may be freely scaled by an arbitrary complex factor. Thus we have a well 



3 In paper I, we regarded u as an unknown rather than solving for it from Eq. (18b); therefore the 
counting was slightly different. 

4 One could alternatively use the algebraic constraint, Eq. (|21]) , for u\(z) to reduce the number of 
complex equations to 5. However, the denominator in Eq. ( f2l] ) vanishes on a region of the negative 
imaginary a axis, and causes numerical trouble. 
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posed eigenvalue problem on a finite domain < z < z 2 which should show a discrete 
spectrum of solutions, apart from linearity. We note that the 1st order system, Eqs. (|19|), 
with its boundary conditions, is not self-adjoint, at least not obviously so; therefore the 
eigenvalue spectrum of a may be expected to be complex. 

Once constructed on the domain < z < Z2, the 1st order solution could be extended 
throughout the spacetime. In the terminology of paper I, Fig. 5, the 1st order solution is 
first constructed in region I, and then it could be extended through regions II and III. We 
will have no need to carry out that extension in this paper. 



V. COORDINATE AND GAUGE CONDITIONS 

Since our field equations possess gauge invariance due to general coordinate invariance 
and global U(l) phase invariance, some unphysical pure gauge modes will appear at 1st 
order, to the extent that the gauge conditions implicit in our boundary conditions Eqs. ( |2"2"| ) 
fail to be unique. 

Indeed there arises a pure gauge mode from an infinitesimal phase rotation — > e le 4> i n 
the Oth order critical solution: 

Si(z) = 0, (27a) 
ui{z) = 0, (27b) 
fa{z)=ifa{z). (27c) 

This gives a time independent solution of Eqs. ( |T§| ) that satisfies the boundary conditions; 
hence it corresponds to an unphysical mode at a = 

A second pure gauge mode arises upon adding an infinitesimal constant to time t — > t + e 
at constant r in the Oth order critical solution. This is possible because our coordinate 
conditions, Eqs. (|5]) normalize t to proper time along the negative time axis (t < 0,z = 0), 
but the zero of time is not specified. Then the solution is perturbed by 

(28a) 
(28b) 

(28c) 

This pure gauge mode has time dependence e~ lar = e~ T and so has negative imaginary 
a = —i. There are no other pure gauge modes in the a plane. These two modes should 
appear as numerical solutions, but are unphysical. 



&i(t, 


*) 


- dbo\ 
dt \ r 


= -(z/t)b'(z) 


= e- T zb'(z), 


ui(t, 


*) 


- - duo 1 
dt \ r 


= -{z/t)u'{z] 


) = e~ T zu'(z), 


01 (T, 




-iujr d(e l " T ip ) 
C dt 


\ r = e~ T (-iuj<j)o[ 


[z) + z(f)' (z)). 



VI. NUMERICAL METHOD 

To solve the 1st order problem we used a Runge-Kutta integrator with adaptive stepsize 
as part of a standard two point shooting method 0, shooting from z = and from z = 
z 2 = 5.0035380, and matching in the middle z = z m = 2.5. For convenience we solved the 
Oth order system, Eqs. flT8|), and the 1st order system, Eqs. ( |I9| ) simultaneously with the 
same steps in z. As discussed in paper I, the similarity horizon z 2 is a demanding place to 
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enforce a boundary condition, and a second order Taylor expansion of the regular solution 
was used for this purpose. 

To solve the 1st order system, we collected all the boundary values but a into a complex 
6-vector X = (s(0), S(0), 61(^2), ^1(^2)) r i z 2)i R( z 2))- Because the equations are linear, the 
matching conditions at likewise linear in the boundary values. A solution is 

found when the values at z m of (pi, Ui, r, s, R, S) upon integrating from z = match with 
those found by integrating from z = Z2, for some boundary values X. We can express this 
matching condition 

A(a)X = (29) 

where A(a) is a 6 x 6 complex matrix which is a nonlinear function of a, constructed 
numerically by integrations of the 1st order equations, Eqs. (fL9D , for six linearly independent 
choices of boundary values X. The condition on a for a solution is then 

det A(a) = 0. (30) 

Once a value for a was found that satisfies this condition, the corresponding boundary values 
X were found as a zero eigenvector of the matrix A; these come in one (complex) parameter 
families, as observed above. Solution of Eqs. (|19|) with boundary values X yields the normal 
mode itself. 

Now, A(a) has been carefully constructed so that it is a complex analytic solution of a. 
This follows from the fact that all equations leading to A contain a but not cr*, together with 
some standard theorems about ODE's. Moreover, A(a) has no singularities in the closed 
lower half a plane. These properties allow us to use a number of ideas from scattering theory 
to study det A(a). In particular, there is a theorem for counting the number Nq of zeros of 
det A within any closed contour C in the closed lower half a plane: 

A c Argdet A = 2rcN c (31) 

where Arg det A is the phase of det A, and A^Arg det A is the total phase wrap (in radians) 
around the closed contour C, a result similar to Levinson's theorem for counting resonances 
in quantum scattering theory. 

Furthermore, there holds a conjugacy relation, 

A*{-a*)=A(a), (32) 

which means that A need only be evaluated for Recx > in the lower half plane. 

The nonlinear equation det A(a) = was solved by the secant variant of Newton's 
method |§. The equation being complex-analytic, the 1-complex-dimensional realization of 
the method was used, and it performed well. 

The following numerical checks are used in our work. First, the two pure gauge modes, 
Eqs. ( p7| , f28|), appeared numerically at the expected places a = 0, —%. Second, the constraint, 
Eq. (PI), was confirmed numerically. Third, we numerically spot-checked the Cauchy- 
Riemann equations for det A(a), to confirm that this function is analytic. Generally the 
accuracy was about 10~ 6 . 
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VII. RESULTS 



There are five values of a in the lower half plane at which det A vanishes, signalling five 
modes. All five are simple zeros; see Fig. 1. 



Recr 

-5 -4 -3 -2 -1 1 2 3 4 5 



1 



oscillatory mode 



jauee mode 



jauge mode 



oscillatory mode 



black hole mode 



4 L Ima 



FIG. 1. The lower half complex a plane. There are five zeros of A(a), representing five modes 
(solutions of the 1st order perturbed field equations). • • Two unphysical gauge modes. ® <g> ® 
Black hole mode, and pair of oscillatory unstable modes. See Sect. VII. 



Two of these modes are pure gauge modes, Eqs. (|27|, ^Bp, which help check internal 
consistency of the work, but are unphysical; as expected these lie on the negative imaginary 
axisQ at a — 0, — i. 

The third solution also lies on the negative imaginary axis, at 

a = -2.58324i ± 0.00002i. (33) 

We interpret this mode as the instability toward forming an infinitesimal^ black hole, Eq. ([[]). 
This corresponds through Eq. (0) to a critical exponent of 7 = 0.387111 ± 0.000003. This 



5 To within Recr = ±1(T 6 . 
6 In classical gravity. 
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value differs from that found by Choptuik jl] for the real scalar field, further evidence of 
non- universality. 

The fourth and fifth solutions form a conjugate pair lying at 

a = ±(3.50224 ± 0.00001) - (0.25867 ± 0.00001)1 (34) 

The presence of this pair of zeros in the lower half plane means that the continuously 
self-similar critical solution is unstable to perturbations, even among solutions "tuned" to 
the black hole threshold. These represent a pair of modes that are not just growing, but 
oscillating as well. These modes must grow until they become nonlinear, and then plausibly 
they go over to Choptuik's solution. 

We evaluated detA(a) around a large rectangular contour C bounded by < Recr < 25, 
> Imcr > — 10 (detouring around the three zeros on the imaginary axis) and used Eq. ( |3~1~D 
to count the zeros lying within. The absence of other zeros in this large rectangle was 
thereby established. Asymptotically, we observed 

detA(a) ~ exp(r oo 0"), 6i (35) 

in both the real and imaginary directions, and the outer boundaries of the large rectangle 
appeared to be well into the asymptotic regime. Thus it appears very likely that there are 
no zeros in the lower half a plane other than the five described above. 

A comment is in order concerning the error of our numerical solution. The error for 
the complex Choptuon constructed in paper I was a little better than 10~ 6 . To investigate 
how the error of the unperturbed solution affected the results in the perturbed problem, we 
varied the parameters of the unperturbed Choptuon by 10 -5 . The effect on the above zeros 
was a variation of about 10~ 5 . We also varied the initial values for the perturbed problem 
(independent of the change in the unperturbed choptuon) in order to study whether the 
locations of the zeros were sensitive to this variation. The resulting change in the values was 
on the order of 10 -6 . This exercise also verified that the Taylor expansions at the endpoints 
and Z2 are of adequate order. 

Since our integrator uses an adaptive stepsize, we did not vary the stepsize directly; 
rather we varied the algorithmic error tolerance used to adjust the stepsize, and we allowed 
the initial stepsize, e, and the location of the matching point, z m , to be random numbers 
within a specific range: 

10 ??? < e < 10~ 6 , 0.5 <z m < 4.5 . (36) 

The combined effect of varying the parameters of the original choptuon by several times 
their quoted error, and of varying the initial values and other parameters for the perturbed 
problem, then gave the overall uncertainties quoted above. 

Thus we conclude that the critical exponent obtained here, 7 = 0.387111 ± 0.000003, is 
indeed different from others that have been calculated to date, such as 7 ra d = 0.3558019, as 
found for radiation fluid collapse by Koike, Hara & Adachi [|J. 

VIII. DISCUSSION AND CONCLUSION 

We have perturbed the continuously self-similar complex choptuon || and find the critical 
exponent for black hole mass, Eq. ([!]) to be 
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7 = 0.387111 ±0.000003, 



(37) 



further evidence for non-universality of 7. We also find a further, oscillatory instability in 
this choptuon. 

In spherically symmetric collapse of a real scalar field, Choptuik JTJ observed that the real, 
discretely self-similar choptuon could be constructed by tuning only one (real) parameter 
p of the infinite number of parameters in the space S of initial conditions. That is, the 
real choptuon is an attractor in a subspace of codimension 1 in <S rea p The results here 
for complex scalar fields imply that the continuously self-similar complex choptuon would 
require tuning of some three (real) parameters p, q, r in the space «S corn pi ex of spherically 
symmetric initial conditions — one (p) to control the real black hole mode, two more (q + ir) 
to control the pair of oscillatory modes. Therefore the complex choptuon appears to be an 
attractor of codimension 3 in 5 com pi ex On threshold — with only p tuned — it is highly 
plausible that the oscillatory instability found here in the complex choptuon would evolve 
nonlinearly toward the original real choptuon. Otherwise, it would have to evolve toward 
still a third choptuon, for which there is currently no evidence. 

Numerical relativity [f[],[9|,|l0"| is needed to follow this nonlinear evolution. In fact, the 
growth rate of the oscillatory mode (measured by — Imcx) is small, ~ 10% of the growth 
rate of the black hole mode, Eqs. ( |33| , |34] ), and consequently the oscillatory instability may 
be hard to isolate numerically. It would be interesting to study related models such as a 
complex scalar field with a conformal coupling £, or the axion/dilaton fields M, to see what 
happens to the oscillatory growing mode. This mode might conceivably move into the upper 
half a plane to become stable. 

Choptuons thus show both similarities and differences compared to stationary black 
holes. For classical black holes, the uniqueness theorems show that "black holes have no 
hair" , and perturbation theory proves that black holes are stable in vacuum, or when coupled 
to massless fields. On the other hand, classical choptuons coupled to a massless complex 
scalar field show at least two alternative states, namely the original choptuon |IJ and the 
complex choptuon ||. In this paper, perturbation theory has shown that the latter state is 
unstable, presumably toward the former. 
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